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Summary 
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presented.  A  modular  approach  is  used  which  permits  each 
tracker  to  optimize  estimation  parameters  for  the  line  assigned. 
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SECTION  I 


Spectral  analysis  procedures  of  various  kinds  are  often  applied  to  what  may  be  loosely 
called  the  "line-tracking"  problem.  If  a  signal  component  of  interest  is  described  as  a  "line", 
this  implies  a  confined  or  "narrow"  spectral  width.  Such  descriptors  are  relative  and  some¬ 
what  vague.  I  remember  well  the  airline  captain  who  reassured  a  nervous  fellow  passenger 

before  takeoff  in  the  bygone  piston/propeller  era  that  all  of  his  flights  were  smooth - but  that 

some  flights  were  smoother  than  others.  In  a  corresponding  fashion  all  of  the  line  signals 
of  interest  here  are  narrow  with  some  lines  being  much  narrower  than  others.  In  spite  of 
the  imprecisionof  definition,  the  line  tracking  operation  may  be  considered  as  one  in  which 
the  parameters  of  signals  qualifying  as  "lines"  are  to  be  determined  and  kept  updated.  Such 
parameters  may  include  center  frequency,  spectral  width  or  shape,  mean  power  level,  and 
relative  phase. 

A  well-known  class  of  trackers  involves  comb  filtering  in  the  frequency  domain  pro¬ 
vided  by  processors  such  as  fast- Fourier  transforms  (FFT’s).  The  individual  comb  signals 
are  detected  and  the  ensemble  of  detected  outputs  is  presented  to  associative  processors 
which  sense  contiguous  amplitude  anomalies  in  the  frequency-time  space  and  declare  a  "track" 
condition  depending  on  program  specifications.  Various  visual  presentations  of  temporal 
spectral-level  (e.g.  lofargrams)  data  have  also  been  used  with  good  success  to  provide 
a  line  detection  and  tracking  capability. 

The  wide  variety  of  line  trackers  which  exist  and  others  which  are  currently  under 
development  give  some  indication  of  the  importance  of  these  techniques  and  of  the  failure  of  a 
"universal"  technique  to  emerge.  Some  of  the  difficulties  of  the  comb  filter,  postdetector 
processor/tracker  approach,  for  example,  are  worthy  of  a  short  discussion.  The  use  of  an 
up-front  coherent  processor  such  as  an  FFT  for  all  lines  requires  that  the  frequency  resolu¬ 
tion  choice  be  a  compromise.  Some  lines  will  span  many  FFT  cells  while  others  may  occupy 
a  very  small  part  of  a  single  frequency  resolution  cell.  In  the  former  case  special  post¬ 
detection  processor  logic  will  be  required  in  order  that  the  signal  smear  (relative  to  the 
spectrum  analyzer)  may  be  recognized  as  a  line.  In  the  latter  case  the  potential  of  the  line 
in  terms  of  detection  level  and  frequency  or  phase  determination  is  squandered  by  an  obvi¬ 
ously  inadequate  processor.  The  compromise  frequency  resolution  choice  for  the  head-end 
FFT  is  seldom  a  satisfactory  one  since  line  widths  encountered  in  may  practical  appli¬ 
cations  can  vary  by  several  orders  of  magnitude. 
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In  some  situations  many  of  the  lines  present  are  of  no  real  interest  except  for  possible 
anomalous  behavior.  Once  these  lines  are  detected  and  identified,  one  only  wishes  to  be 
made  aware  of  an  unusual  event  such  as  the  sudden  disappearance  of  one  or  more  of  these 
lines.  During  routine  periods  such  lines  act  really  as  jammers,  cluttering  up  displays  and 
interfering  with  signal  processing  operations  on  other,  often  much  weaker  signals  which  may 
be  of  primary  interest.  The  frequency-quantization  property  of  comb  filters  such  as  the 
FFT  is  one  of  the  prices  paid  for  economy  in  spectral  computation.  In  the  tracking  operation 
such  compartmentization  by  frequency  can  be  troublesome  since  a  line  signal  will  often  move 
through  several  frequency  cells  in  the  course  of  an  analysis.  Computation  of  line  frequency 
and  relative  phase  can  also  be  complicated  by  the  combination  of  a  signal  which  exists  in  a 
frequency  continuum  and  an  analysis  procedure  which  imposes  the  nuisance  of  frequency 
quantization. 

The  nature  of  the  line  tracking  problem  varies  greatly  with  specific  application.  The 
residual  signal  analysis  (RSA)  procedure  presented  here  is  a  technique  which  has  some  unique 
and  very  useful  properties.  In  some  applications  the  RSA  subsystem  can  satisfy  the  entire 
line  tracking  requirement.  In  other  cases  these  algorithms  can  work  best  as  a  preprocessor 
ahead  of  existing  spectral  analysis  systems. 

1.  BASIC  PRINCIPLES 

Operation  of  the  residual  signal  analysis  (RSA)  system  will  now  be  described  with  ref¬ 
erence  to  Figure  1-1.  Both  time  and  level  quantization  are  assumed  in  the  descriptions 
which  follow.  Three  main  signal  lines  or  buses  are  used:  the  Input  Signal  Bus(ISB),  the 
Tracker  Signal  Bus(TSB)  and  the  Residual  Signal  Bus(RSB).  Discrete  processing  entities 
called  tracker  modules  (TMf)  are  made  available  and  are  assigned  exclusively  to  each  line 
being  tracked.  These  processing  entities  may  be  instrumented  in  hardware,  software  or  in 
combination.  The  methodology  internal  to  the  tracker  modules  can  be  of  various  types  which 
will  be  described  and  discussed  later. 

For  the  moment  it  is  sufficient  to  understand  that  a  tracker  module  is  assigned  to 
one  spectral  line  only.  The  module  contains  logic  which  senses  that  the  assigned  line  has 
indeed  been  "acquired"  so  that  an  internal  "track"  flag  may  be  set  true.  An  external  switch 
Sr  is  closed  only  when  "track"  is  true,  otherwise  the  switch  is  open  and  that  module  makes 
no  contribution  to  the  Tracker  Signal  Bus 

Once  track  is  achieved  the  tracker  module  optimizes  its  operating  parameters  to  match 
the  spectral  width  and  dynamics  of  the  assigned  line.  After  each  clock  cycle  the  tracker 
module  provides  at  its  output  an  estimate  of  the  signal  value  which  will  exist  at  the  next 
clock  or  sampling  time.  This  predicted  signal  voltage  estimate  is  delayed  (saved)  for  one 
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clock  or  sampling  time.  This  predicted  signal  voltage  estimate  is  saved  so  that  it  will  be 
immediately  available  at  the  next  clock  cycle. 
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Figure  1-1.  Residual  Signal  Analysis  (RSA)  System  Block  Diagram 

At  the  next  sampling  time  the  Input  Signal  Bus  voltage  has  subtracted  from  it  the  Tracker 
Signal  Bus  voltage  which  contains  the  sum  of  signal  estimates  of  operative  trackers.  The 
difference  voltage  is  placed  on  the  Residual  Signal  Bus,  which,  if  all  has  gone  well,  contains 
greatly  attenuated  a  versions  of  the  signals  being  tracked.  This  Residual  Signal  Bus  voltage 
feeds  the  postprocessing  subsystems  (if  any)  and  also  feeds  each  of  the  trackers. 

At  this  point  we  have  a  bit  of  a  problem.  If  the  trackers  are  doing  a  good  job  then  there 
will  be  only  vestiges  of  the  input  spectral  line  signals  on  the  Residual  Signal  Bus  which  forms 
the  common  drive  to  the  trackers.  There  will,  in  fact,  be  very  little  in  the  way  of  lines  for 
the  trackers  to  track.  If  track  is  lost,  the  switches  Sr  open,  the  Input  Signal  Bus  voltage  and 
the  Residual  Signal  Bus  voltage  become  identical  and  the  spectral  lines  reappear  at  the 
tracker  inputs. 
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In  order  to  resolve  this  difficulty,  feedback  is  provided  around  each  tracker  module. 
This  at  first  appears  to  be  positive  feedback  but  in  fact,  it  is  not.  Note  in  Figure  1-1  that 
tracker  module  1,  for  example,  contributes  a  voltage  6 j  to  the  Tracker  Signal  Bus(track  con¬ 
dition  assumed,  switch  Sj  closed).  This  voltage  component  when  subtracted  from  the  Input 
Signal  Bus  voltage  eliminates  the  spectral  line  being  tracked  by  module  1  from  appearing  on 
the  Residual  Signal  Bus.  If  this  voltage  6^,  is  added  back  to  the  Residual  Signal  Bus  voltage 
to  form  the  actual  input  or  drive  to  tracker  module  1,  as  shown,  then  there  is  in  fact  no 
feedback  around  the  module.  Note  also  that  the  spectral  line  being  tracked  by  module  1 
appears  unaltered  at  the  input  to  that  module  but  does  not  appear  at  any  other  module  input 
because  of  its  (relative)  absence  from  the  Residual  Signal  Bus. 

If  this  argument  is  extended  to  include  simultaneous  operation  of  a  number  of  trackers, 
it  may  be  seen  that  each  tracker  sees  its  own  signal  unattenuated.  However,  each  tracked 
signal  is  (nearly)  transparent  to  all  of  the  other  trackers  or  other  processors  (such  as  dis¬ 
plays)  that  may  be  fed  from  the  Residual  Signal  Bus.  The  system  is  thus  ’’failsafe"  in  the 
sense  that  a  signal  line  lost  or  dropped  from  track  by  an  assigned  module,  appears  at  once 
on  the  Residual  Signal  Bus.  Generally  speaking  lines  not  seen  on  the  Residual  Signal  Bus 
are  lines  already  accounted  for  by  the  tracker  module  processors.  The  Residual  Signal 
Bus  contains  at  all  times  those  signals  which  remain  after  tracker-module  processing  and 
are  "residual"  to  this  process. 

Physical  signal  elimination  has  many  advantages,  some  of  which  are  not  obvious  at 
first  glance.  Consider  the  original  line  acquisition  problem  relative  to  Figure  1-1.  The 
tracker  control  unit  (TCU)  dispatches  the  first  tracker  (Say  TM^)  to  search  for  and  acquire 
a  spectral  line  within  a  specified  passband.  Once  this  tracker  indicates  acquisition  to  the 
TCU  a  second  tracker  (Say  TM2)  may  be  dispatched.  Since  the  line  already  being  tracked  by 
TMj  is  transparent  to  the  searching  TMg,  this  second  tracker  will  usually  pass  by  the  spec¬ 
tral  line  already  acquired,  and  lock  on  to  a  spectral  line  beyond  this  point.  Upon  receipt  of 
a  lock  signal  from  TM2,  TMg  is  dispatched,  and  so  forth.  The  "search  and  destroy"  nature 
of  this  procedure  can  thus  be  used  to  good  advantage  in  the  initialization  process. 

A  spectral  line  parameter  of  routine  interest  is  line  "center"  frequency.  Such  fre¬ 
quency  estimates  can  be  badly  biased  by  the  presence  of  strong  lines  in  the  frequency 
neighborhood  of  the  spectral  line  of  real  interest.  While  measurement  means  can  be  devised 
which  attenuate  off-frequency  signals  over  extended  skirt  frequency  regions,  such  approaches 
tend  to  be  expensive.  In  the  RSA  system  a  line  tracked  is  a  line  eliminated  and  this  normally 
represents  an  inexpensive  means  for  getting  (effectively)  skirt  attenuation  at  precisely  those 
frequency  locations  where  attenuation  is  needed.  Generally  speaking,  a  strong  line  should 
be  tracked  even  if  it  is  of  no  interest  since  it  can  interfere  with  the  processing  of  other  lines. 
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The  reduction  of  dynamic  range  by  strong  line  elimination  through  RSA  acquisition  can 
be  very  useful.  This  is  especially  true  when  limited  amplitude  tolerance  is  encountered  in  a 
subsystem  such  as  an  amplitude-frequency-time  display.  Much  new  information  can  often  be 
seen  on  such  displays  in  the  region  of  a  strong  line,  once  that  line  has  been  RSA-removed. 

In  fact  the  man-machine  process  can  be  iterative.  Lines  observed  are  removed  by  operator 
module  assignments.  Once  this  is  done,  display  levels  may  be  readjusted  to  reveal  new 
lines  which  in  turn  may  be  given  special  RSA  attention.  The  property  that  the  operator  sees 
only  what  the  hardware  is  not  accounting  for  provides  for  some  very  useful  man-machine 
team  characteristics. 

The  above  qualitative  discussion  has  presented  the  basic  principles  of  the  RSA  tech¬ 
nique  and  has  allowed  some  of  the  advantages  of  this  approach  to  be  indicated.  In  the 
following  sections  some  simple  quantitative  results  are  obtained  and  some  comparisons 
of  tracker  module  implementation  schemes  are  given. 


SECTION  II 

ANALYTICAL  DETAILS 

A  simple  analysis  of  tracker  module  operation  will  now  be  given  with  the  aid  of 

Figure  2-1.  It  is  assumed  that  there  are  K-operative  modules  and  only  the  r1*1  module 

is  shown  in  Figure  2-1.  A  linear  operator  H  (Z)is  assumed  as  the  processor  of  the  module 

r  ij. 

drive  signal  d  to  produce  the  signal  6  which  is  an  estimate  of  the  r  spectral  line 
r,  n  r,n 

voltage  as  it  will  appear  at  the  next  sampling  time  index  (n  +  1).  Specific  implementations 
of  the  module  linear  operator  Hr(Z)will  be  discussed  later.  For  the  present  the  general 
representation  indicated  in  Figure  2-1  is  appropriate. 


Figure  2-1.  Tracker  Module  Details 


The  algebra  of  Figure  2-1  follows  easily.  The  residual  signal  is  the  difference  be¬ 
tween  input  signal  and  Tracker  Signal  Bus  value 


2-1 


t 


Equation  (2-6)  yields 
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where  G(Z),  as  defined  here,  becomes  the  transfer  function  of  the  RSA  subsystem.  Equations 
(2-7)  and  (2-9)  combine  to  give 
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which  when  substituted  into  Equation  (2-8)  yields 
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which  may  be  combined  with  Equation  (2-10)  to  give 
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In  the  special  case  of  a  single  operative  module  (which  will  be  examined  in  what  immediately 
follows)  Equation  (2-13)  becomes 


G(Z)  =  1  -  Z_1  Hj(Z).  (2-14) 

An  estimator/predictor  structure  that  has  given  excellent  results  is  sketched  in 
Figure  2-2.  This  modified  bandpass  filter  may  be  used  as  the  Hr(Z) linear  operator  of 
Figure  2-1.  Detections  along  the  I  and  Q  axes  are  performed  at  the  filter  radian  center 
frequency  and  the  resulting  quadrature  components  are  identically  low-pass  filtered. 
These  recovered  i  and  Q  components  remodulate  carrier  signals  that  are  advanced  in  phase 
by  uqT,  the  carrier  phase  shift  in  one  sampling  period.  The  output  signal  6(t)  is  an  estimate 
of  the  spectral  line  signal  as  it  will  appear  T-seconds  in  the  future.  Advantage  is  taken  of 
the  fact  that  the  narrowband  I  and  Q  components  of  the  spectral  line  will  change  negligibly 
during  the  sampling  interval  T  and  that  only  the  carrier  phase  need  be  advanced  to  achieve 
good  prediction.  Since  the  sampling  rates  in  this  system  must  accomodate  the  total  analysis 
bandwidth  and  since  the  spectral  line  widths  are  very  small  compared  to  the  analysis  band¬ 
width,  the  approximation  involved  here  is  a  very  good  one. 
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(R,  C,  OR  L  OF 
FIGURE  3-2 ) 

Figure  2-2.  Estimator/Predictor  Arrangement  for  Narrowband  Signals 

This  linear  operator  is  further  simplified  by  specifying  that  the  low-pass  filters  are 
single-pole  filters  of  variable  bandwidth.  The  input  (xR)  and  output  (yR)  samples  of  the 
low-pass  filters  of  Figure  2-2  satisfy 

yR  =  (l-o)  x„  +«yn_1  (2-15) 

where  the  3-dB  cutoff  frequency  fu  is  related  to  the  parameter  cl  by 

ll 

-2^f  T 

a  =  e  "  (2-16) 

The  3-dB  bandpass  bandwidth  of  the  Figure  2-2  filter  is  as  a  consequence  2f„. 

The  use  of  a  simple  single-pole  filter  in  place  of  a  more  complicated  structure  may  be 
defended  for  most  applications  without  difficulty.  The  detailed  spectral  shape  of  the  line 
signals  is  seldom  known  precisely  so  that  "optimum"  processing  is  not  possible.  Even  if  the 
power  density  spectrum  were  known  with  precision,  a  single-pole  filter  of  appropriate  band¬ 
width  will  usually  perform  so  close  to  "optimum"  that  the  difference  is  swamped  by  the  errors 
inherent  in  the  many  other  modeling  assumptions  made. 

It  is  convenient  at  this  point  to  consider  "complex  envelope"  representations  of  the 
narrowband  signals.  The  mathematical  development  of  Equations  (2-1)  through  (2-14)  can 
apply  to  both  real-  and  complex-envelope  representations  of  signals  if  appropriate  care  is 
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taken.  If  the  frequency  term  implied  by  Z  =  exp(j2irfT)  is  now  measured  from  reference 
frequency  ^(cu Q  =  2n{^)t  the  transfer  function  of  the  processor  of  Figure  2-2  becomes 


(1  -  a)Z 
(Z  -cO 


(2-17) 


The  additional  delay  of  T  seconds  implied  by  the  Z  operator  of  Figure  2-2  requires  that 
the  Z-transform  of  the  complex  envelope  of  the  undelayed  signal  be  multiplied  by  Z-1  and 

-KT 

e  to  account  for  carrier  phase  rotation.  When  this  is  done  for  Equation  (2-17)  and  the 
result  used  in  place  of  the  Z  1H1(Z)  term  of  Equation  (2-14),  one  obtains  the  RSA  subsystem 
transfer  function  for  a  single-pole  estimator/predictor  based  on  a  complex  envelope  signal 
representation.  Using  Equations  (2-14)  and  (2-17)  in  this  manner  yields 
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The  phase  response  from  the  input  signal  bus  to  the  residual  signal  bus  is  shown  in 
Figure  2-3.  The  phase  is  odd-symmetric  about  f  =  0,  hence  there  is  a  180°  jump  in  phase 
at  f  =  0.  The  frequency  response  curve  from  input  signal  buB  to  residual  signal  bus  is  also 
shown  in  Figure  2-3.  A  symmetric  notch  results  with  infinite  attenuation  at  band  center  (f  =  0). 
This  implies  that  a  pure  CW  signal  would  be  perfectly  predicted  and  thus  entirely  removed 
from  the  residual  signal  line.  This  represents  neither  much  of  a  surprise  nor  much  of  an 
accomplishment.  The  behavior  of  this  configuration  for  more  realistic  signal  types  will  be 
examined  in  the  next  section.  Note  from  Figure  2-3  that  at  the  3  dB  frequencies  of  the  single¬ 
pole  filter  an  insertion  loss  of  3-dB  in  the  G-function  exists.  This  gives  some  indication  of 
the  "spillover"  effects  in  the  frequency  domain  of  a  tracker  module  of  this  type.  At  these 
frequencies  the  phase  shift  will  be  45°. 

The  estimator/predictor  operator  discussed  above  forms  one  part  of  a  particular 
tracker  module  configuration.  The  complete  tracker  module  arrangement  will  be  discussed 
in  the  next  section  along  with  a  consideration  of  more  realistic  signal  spectra. 
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SECTION  III 


TRACKER  MODULE  CONFIGURATIONS 

The  phase-lock  loop  arrangement  of  Figure  3-1  has  been  used  with  success  as  a  tracker 
module  in  the  context  of  Figure  1-1.  In  spite  of  the  fact  that  both  I  and  Q  detections  are 
made,  this  is  not  a  Costas  loop.  It  is  a  normal  carrier-lock  loop  with  the  Q  channel  providing 
the  error  signal.  The  signal  envelope  is  obtained  from  the  I  and  Q  signals  and  is  used  as  a 
threshold  against  which  the  average  component  of  I  is  tested.  If  the  smoothed  I  signal  is 
sufficiently  close  to  the  smoothed  envelope  value  "lock"  is  declared  and  the  logic  line  goes 
true"  and  enables  or  closes  switches  S.  The  voltage  controlled  oscillator  (VCO)  signals  are 
advanced  in  phase  by  w0T  and  are  remodulated  by  I  and  Q  as  shown  to  obtain  the  predicted 
output  signal  0(t).  The  similarities  between  Figures  3-1  and  2-2  are  clearly  evident.  How¬ 
ever  the  operational  philosophies  from  which  these  two  figures  are  drawn  are  quite  different. 
Figure  3-1  implies  an  actual  phase  lock  to  the  signal  with  the  temporal  wanderings  of  signal 
phase  being  continually  determined  and  followed.  Figure  2-2  implies  a  filtering  arrangement 
which  expects  that  center  frequency  and  bandwidth  will  be  sensibly  selected,  but  no  demand 
is  made  with  regard  to  phase  coherence. 


Figure  3-1.  Phase-Lock  Loop  Tracker  Module  Arrangement 
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Figure  3-2  shows  a  three  filter,  frequency  discriminator  loop  tracker  module  that  has 
given  very  good  results.  The  filter  structures  are  of  the  type  shown  in  Figure  2-2.  The 
square-law  detected  outputs  of  the  outboard  filters  are  used  to  obtain  a  frequency  tracking 
error  signal  through  the  computation  of  (R  -  L)/R  +  L).  This  error  signal  enters  a  filter 
frequency  control  unit  which  combines  the  direct  error  signal  with  its  integral  to  produce  a 
filter  frequency  (group)  control  signal,  which  is  initialized  by  some  means  to  represent  the 
desired  spectral  line  radian  frequency  w  .  Offset  control  signals  are  also  provided  for  the 
outboard  filters  as  shown.  In  the  examples  of  this  paper  the  filter  center  frequencies  are 
spaced  by  2fH>  the  3-dB  single-pole  passband  bandwidth.  When  a  frequency  error  occurs, 
the  error  signal  causes  the  three  filters  to  move  together  such  that  the  center  filter  tracks 
the  spectral  line. 
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Figure  3-2.  Frequency  Discriminator  Loop  Tracker  Module  Arrangement 
Track  indication  is  derived  from  the  Q-function  calculation 

Q  =  (3-1) 

which  displays  an  interesting  and  useful  variation  with  filter  bandwidth  (kept  the  same  for 
all  three  filters)  for  a  given  signal  spectrum  and  background  noise  level.  If  the  filter  band- 
widths  are  made  very  large  compared  to  signal  bandwidth,  the  output  power  levels  of  the 
filters  (which  is  what  the  smoothed  R,  C,  and  L  represent)  approach  equality.  This  occurs 
because  background  noise  power  tends  eventually  to  dominate  signal  power  as  the  filter 
bandwidths  are  increased.  When  this  occurs  the  Q  function  of  Equation  (3-1)  approaches 
zero. 
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At  the  other  extreme  consider  filter  bandwidths  that  are  very  narrow  compared  to  sig¬ 
nal  bandwidth.  The  combined  signal  and  noise  spectrum  is  now  essentially  constant  and  equal 
in  the  response  regions  of  the  three  filters  so  that  the  R,  C,  and  L  levels  again  become  equal. 
This  causes  the  Q  function  to  tend  once  again  toward  zero.  The  Q  function  can  be  expected 
to  have  a  peak  value  for  some  "optimum"  filter  bandwidth.  The  Q  function  is  tested  in  the 
arrangement  of  Figure  3-2  in  order  to  determine  whether  the  "track"  is  true  of  false.  If  it 
is  decided  that  a  line  is  being  tracked,  switch  S  is  closed  and  the  predicted  signal  from  the 
center  filter  (see  Figure  2-2)  is  made  available  on  an  output  line. 

Both  the  phase-lock  loop  (PLL)  tracker  module  of  Figure  3-1  and  the  frequency  dis¬ 
criminator  loop  (FDL)  tracker  module  of  Figure  3-2  have  been  used  satisfactorily.  The 
choice  between  these  two-tracker  module  types  will  depend  on  many  factors  including  the 
specific  application  at  hand.  The  FDL  tracker  is  chosen  for  further  analysis  in  this  paper 
for  several  reasons.  Very  weak  lines  can  be  tracked  using  the  FDL  by  long-term  smoothing 
of  R,  C,  and  L  which  cannot  be  handled  at  all  by  the  PLL  tracker.  This  is  because  the  tem¬ 
poral  behavior  of  the  I  and  Q  components  are  completely  masked  by  noise.  Even  at  higher 
signal  levels  the  signal  envelope  fluctuations  and  the  attendant  phase  discontinuities  that 
accompany  a  fade,  cause  annoying  and  recurring  losses  of  lock  in  the  PLL  case.  Since 
acquisition  is  always  more  difficult  than  tracking,  the  repeated  supervisory  demands  on  the 
Tracker  Control  Unit  (Figure  1-1)  because  of  PLL  dropouts  becomes  a  burden.  In  certain 
signal-to-noise  ratio  (SNR)  regimes  the  PLL  may  give  somewhat  better  estimates  of  signal 
frequency  and  phase,  however  these  quantities  are  certainly  available  from  the  FDL.  Note 
in  Figure  2-2  that  an  arctangent  operation  on  the  smoothed  I  and  Q  components  will  yield 
signal  phase.  It  should  also  be  remembered  that  a  linear  bandpass  filter  with  an  intelligently- 
selected  bandwidth  provides  spectral  line  signal  estimation  performance  which  is,  by  any 
practical  measure,  at  the  "optimum"  performance  level. 

With  regard  to  the  FDL  tracker  module  a  means  for  selecting  bandwidth  will  now  be 
considered.  For  RSA  purposes  estimation  and  prediction  performance  are  essentially  one 
and  the  same  due  to  the  high  ratio  of  sampling  rate  to  signal  bandwidth  as  previously 
discussed. 

For  the  single-pole  filter  of  Equation  (2-15)  if  Z-transforms  are  taken  and  the  usual 
narrowband  approximations  are  made,  the  filter  transfer  function  may  be  represented  by 


where  f  is  measured  from  band  center  and  is  one-half  the  total  bandpass  3-dB  bandwidth. 
For  convenience  in  discussion  let  the  spectral  line  signal  have  a  (two-sided)  power  density 
spectrum  of  unity  over  a  bandwidth  of  B  Hz.  The  signal  to  noise  ratio  will  be  indicated  by  y . 

If  the  mean-square  error  4>  is  taken  as  the  measure  of  the  goodness  of  signal  estimation, 
then  by  the  usual  calculus  one  obtains: 


(3-3) 


The  results  are  shown  in  Figure  3-3.  The  filter  bandwidth  required  for  best  estimation  varies 
with  SNR.  Two  types  of  errors  exist:  noise  and  signal  distortion.  Increasing  bandwidth  in¬ 
creases  noise  but  lowers  signal  distortion.  The  optimum  bandwidth  increases  with  increasing 
SNR. 

It  will  now  be  shown  that  the  filter  bandwidth  which  maximizes  the  Q-function  (a  readily 
measurable  quantity)  will  also  result  in  near-optimum  estimation  performance.  As  before, 
the  signal  spectrum  is  assumed  to  have  a  (two-sided)  unit  power  density  spectrum  over  a 
passband  bandwidth  of  B  Hz.  The  three  single-pole  filters  have  a  bandwidth  of  2fR  and  a  center 
frequency  spacing  of  2fR.  Letting  x  =  fH/B,  V  =  SNR,  routine  analysis  yields 
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The  results  are  shown  plotted  in  Figure  3-4.  Comparing  Figures  3-3  and  3-4  shows  that  a 
bandwidth  adjustment  for  maximum  Q  yields  nearly  optimum  estimation  performance.  While 
this  simple  analysis  has  idealized  the  signal  spectrum  and  used  mean-square  error  as  the 
performance  measure,  practical  experience  has  indicated  that  adapting  to  peak  Q-values  is 
indeed  a  viable  strategy. 


As  a  practical  matter  parameter  optimization  has  not  proved  to  be  a  critical  problem. 
The  performance  curves  in  the  real  world  appear  to  have  broad,  smooth  peaks.  Furthermore, 
the  nonstationary  nature  of  most  signals  and  noises  encountered  makes  precise  performance 
peaks  very  difficult  to  find  and  verify  in  any  event. 
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SECTION  IV 


SEA  TEST  EXAMPLES 


The  data  processed  in  the  examples  which  follow  was  obtained  from  completely  unclass¬ 
ified  sources.  The  test  scenario  involved  transmission  from  one  vessel  of  four  highly  stable 
and  coherent  spectral  lines  at  187  Hz,  195  Hz,  200  Hz,  and  206  Hz  (’'projector  lines").  A 
broadband  noise  signal  was  also  radiated  by  the  transmit  vessel.  These  signals  were  re¬ 
corded  aboard  another  vessel.  When  the  recordings  were  analyzed  a  fifth  line  was  discovered 
near  215  Hz.  This  line  was  found  Lo  be  generated  by  a  pump  aboard  the  receive  vessel  and  is 
described  in  what  follows  as  a  "pump  line". 

The  transit  of  one  vessel  past  the  other  during  the  test  caused  the  four  projector  lines 
to  display  a  doppler-time  history  pattern  of  interest.  The  task  of  the  RSA  processor  sub¬ 
system  involved  accurate  measurement  of  this  Doppler-time  history  for  the  four  projector 
lines.  In  addition  it  was  required  that  all  lines  (including  the  pump  line)  be  removed  to  pro¬ 
vide  access  to  the  broadband  noise  signal  which  was  also  radiated  by  the  transmit  vessel. 

This  broadband  signal,  unencumbered  by  narrowband  line  "interference",  was  to  be  delivered 
to  a  post-processing  subsystem  connected  to  the  Residual  Signal  Bus  (Figure  1-1).  In  the 
following  the  frequency  measurement  and  then  the  line-removal  performance  of  the  RSA  sub¬ 
system  will  be  discussed. 

In  Figure  4-1  plots  are  shown  of  the  measured  projector  line  frequencies.  This  data 
is  taken  from  computer  listings  which  contained  the  measured  frequencies  listed  to  the  near¬ 
est  millihertz.  Measurements  were  printed  every  5s  and  the  time  span  of  the  data  of  Figure 
4-1  was  chosen  to  include  the  closest  point  of  approach  (CPA)  of  the  two  vessels  which  occur¬ 
red  at  a  range  of  about  4000  yards.  The  filter  bandwidths  for  the  four  modules  assigned  to 
the  projector  lines  were  100  mHz.  A  fifth  module  was  assigned  to  the  pump  line  at  215  Hz 
with  a  filter  bandwidth  of  1  Hz.  This  last  tracker  was  assigned  for  line-elimination  purposes 
only,  the  frequency  data  from  this  module  was  ignored. 

An  estimate  of  the  quality  of  the  frequency  measurements  of  Figure  4-1  was  obtained 
by  the  process  to  be  described.  If  the  transmit  frequency  for  line  n  is  represented  by  F^-n 
and  the  received  frequency  by  FRn,  then 
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where  v  is  the  relative  closing  velocity  between  the  two  vessels  and  c  is  the  sound  velocity 
in  the  medium.  One  line  may  be  used  as  a  reference  (line  2  in  this  case)  to  obtain  a  (1  +  v/c) 
value  using  Equation  (4-1).  The  other  line  frequencies  may  then  be  predicted  from  this  value 
using  the  known  transmit  frequencies.  Let  En  represent  the  error  between  the  predicted 
frequency  and  the  measured  frequency.  Then 


E  =  — —  F* 
n  Ft2  r  R2 


F’Rn,  n  =  l,  3.4 


where  F'R2  and  F'Rn  are  the  measured  received  frequencies. 


(4-2) 


If  the  received  frequency  measurements  are  assumed  to  be  unbiased,  these  measure¬ 
ments  may  be  represented  by 


F'b„  -  FR„  ♦ 


(4-3) 


where  FRr  represents  the  true  signal  frequency  at  the  receiver  and  6  represents  a  zero- 
mean  error  term  in  the  measurement.  Equation  (4-2)  may  be  rewritten  as 


f  =  _Ia 
n  "  ft2 

However,  since 


(FR2  +  6 2)  ~  (FRn  +  6n)  *  n  ~  3.  4  . 


(4-4) 


R2 
'  T2 


F  =  (l  +  —  ^  F™  =  F  ,  n  =  1,  3,  4 

rTn  \  o  I  rTn  rRn  ’  *  ’ 


(4-5) 


then 


Jn  V  FT2  ) 


&2  -  6n  ,  n  =  1,  3,  4 


(4-6) 


The  error  quantities  En  should  therefore  have  zero  mean.  As  measured  for  the  run  of 
Figure  4-1,  there  resulted  Ej  =  -0.05,  Eg  =  +0.18,  and  =-0.30  mHz,  giving  an  over¬ 
all  average  of  -0.06  mHz.  At  time  955  in  this  run  the  measured  frequencies  as  listed  were 
exactly  187.000,  195.000,  200.000  and  206,000  Hz.  These  are  precisely  the  synthesized 
transmit  frequencies.  Since  this  time  lies  within  the  exercise  track  reconstruction  CPA 
window,  we  have  further  evidence  that  the  frequency  measurements  produced  by  the  RSA 
modules  are  unbiased. 
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If  the  measurement  error  terms  6  are  assumed  to  be  independent,  zero  mean  and 

2 

have  the  same  mean-squared  value  cr^  ,  squaring  both  sides  of  Equation  (4-6)  and  averaging 
yields 

n  =  1,  3,  4  (4-7) 


The  three  values  of  ^En  J  of  the  Figure  4-1  data  were  0.66,  1.46,  and  1. 21  which  yields 

values  of  0.34,  0.71,  and  0.57.  Averaging  these  results  gives  a  c r  of  0.73  mHz. 

6  o 

These  results  indicate  that  the  frequency  measurements  for  the  run  of  Figure  4-1  are 
essentially  unbiased  and  have  a  standard  deviation  of  the  order  of  three-quarters  of  one 
millihertz. 

Figure  4-2  shows  a  spectrogram  of  the  RSA  output  (Residual  Signal  Bus)  signal  with 
one  operative  tracker  module  assigned  to  the  pump  line  at  215  Hz.  As  a  consequence  this 
line  does  not  appear  in  Figure  4-2,  however  the  four  projector  lines  are  clearly  visible  at 
187,  195,  200  and  206  Hz.  The  broadband  noise  signal  radiated  by  the  transmit  vessel 
extends  from  about  100  Hz  to  500  Hz  and  the  noise  background  peak  in  the  vicinity  of  150  Hz 
is  believed  to  arise  from  hull  noise  at  the  receive  vessel. 

Figure  4-3  shows  a  spectrogram  of  the  RSA  signal  output  when  five  trackers  are  de¬ 
ployed,  one  for  each  of  the  projector  lines  and  one  for  the  pump  line.  (Plot  normalization 
results  in  a  +5.4  dB  shift  of  Figure  4-3  data  relative  to  that  of  Figure  4-2. )  As  may  be 
seen,  all  five  dominant  lines  have  been  suppressed  so  that  background  noise  and  other 
narrowband  components  now  dominate  the  spectral  plot.  In  these  last  two  figures  an  (FFT) 
window  of  one  second  was  used  and  five  FFT  outputs  were  averaged  to  produce  the  data 
shown.  The  time  periods  from  which  the  data  were  drawn  were  about  2. 5  minutes  apart. 

A  different  perspective  on  frequency-domain  effects  may  be  obtained  from  the  lofar- 
grams  shown  in  Figures  4-4  to  4-6.  The  first  of  these  amplitude  versus  frequency  time- 
history  displays,  Figure  4-4,  shows  the  RSA  input  signal.  The  five  lines  dominate  the  dis¬ 
play  as  expected.  Figure  4-5  shows  the  RSA  output  signal  when  five  trackers  are  deployed. 
The  background  noise  peak  in  the  150  Hz  neighborhood  now  saturates  the  display  as  the  five 
line  levels  have  been  significantly  reduced.  Vestiges  of  the  four  projector  lines  may  still 
be  seen  in  Figure  4-5  even  through  their  level  appears  to  be  reduced  below  that  of  the  back¬ 
ground  noise  peak  near  150  Hz.  Note  that  the  pump  line  (compare  with  Figure  4-4)  has  been 
reduced  to  a  broad  smear  in  Figure  4-5,  This  pump  line  is  much  more  typical  of  real-world 
lines  than  the  projector  lines,  of  course. 
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Figure  4-5.  Lofargram  of  RSA  Output  Signal  with  Five  Trackers  Deployed 
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(The  strength,  source  purity  and  short  water  path  of  the  projector  lines  create  some 
unusual  problems  in  the  demonstration  of  the  RSA  technique.  Estimation  module  bandwidths 
of  100  mHz  were  used  for  the  frequency  measurements  of  Figure  4-1.  This  yields  excellent 
frequency  estimation  but  poor  line  elimination  performance.  The  spectrograms  of  Figures 
4-2  and  4-3  and  the  lofargrams  of  Figures  4-5  and  4-6  were  obtained  using  1  Hz  filter  band- 
widths.  This  gave  much  better  suppression  but  higher  variance  frequency  data.  Wider 
filter  bandwidths  could  have  been  used  to  obtain  much  more  impressive  line  suppression 
performance  for  the  lofargrams.  It  is  strongly  indicated  that  energy  arrival  at  the  receive 
vessel  contains  two  components.  A  carrier  component  is  present  whose  average  frequency 
is  truly  responsive  to  radial  platform  velocity.  Water  path  phase  and  amplitude  modulation 
effects  create  sideband  energy  as  a  second  received  component  which  is  at  a  non-trivial 
level  relative  to  background  noise  due  to  the  very  high  projector  line  levels.  This  creates 
the  unusual  condition  in  which  poor  suppression  accompanies  good  frequency  estimation. 

This  rather  artificial  situation  may  have  much  more  scientific  value  as  regards  medium- 
effects  studies  using  RSA  than  engineering  utility. ) 

Of  special  interest  in  Figure  4-5  is  the  new  line  now  visible  between  projector  lines  3 
and  4.  An  enlargement  of  this  region  of  Figure  4-5  is  shown  in  Figure  4-6.  The  dynamic 
range  reduction  provided  by  the  RSA  subsystem  allows  much  more  low-level  data  to  be  in¬ 
vestigated  using  limited  dynamic  range  displays  as  evidenced  by  Figure  4-6.  While  these 
projector  lines  are  atypical,  the  enhanced  capability  of  the  lofargram  display  when  operated 
in  conjunction  with  an  RSA  subsystem  is  a  routine  and  very  useful  benefit.  The  "new"  line 
uncovered  in  Figures  4-5  and  4-6  stays  between  projector  lines  3  and  4  for  the  entire  run. 

The  etiology  of  this  spectral  component  is  not  understood  at  this  time. 

The  sea  data  processed  here  was  readily  available  and  its  use  was  convenient  for  dem¬ 
onstration  of  the  RSA  technique.  The  highly  accurate  measurements  of  frequency  obtained 
were  made  possible  by  the  unusual  quality  of  the  source  of  the  projector  lines  and  by  the 
short  water  path  involved  in  acoustic  propagation.  In  most  real-world  situations,  conditions 
are  not  as  fortuitous  as  those  found  here,  of  course,  and  somewhat  more  uncertainty  in  fre¬ 
quency  estimation  is  inevitable.  Since  individual  tracker  module  parameters  may  be  tailored 
to  the  specific  characteristic  of  the  assigned  line,  frequency  measurement  may  be  optimized 
for  each  line  of  interest.  The  absolute  level  of  optimized  frequency  estimation  will  vary  con¬ 
siderable  from  line  to  line  depending  upon  spectral  purity,  signal-to-noise  ratio,  and  center 
frequency  temporal  behavior  (line  dynamics). 
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SECTION  V 


CONCLUSIONS 

The  residual  signal  analysis  technique  represents  a  unique  and  unusual  approach  to 
the  line  tracking  problem.  Because  of  its  novelty  and  the  limited  resources  committed  thus 
far  to  its  development,  much  more  work  remains  to  be  done  in  order  to  exploit  the  inherent 
advantages  of  the  concept. 
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